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<N ■ ABSTRACT 

^ ■ We highlight the potential importance of gaseous TiO and VO opacity on the highly irradiated close-in giant 

Q planets. The atmospheres of these planets naturally fall into two classes that are somewhat analogous to the M- 
and L-type dwarfs. Those that are warm enough to have appreciable opacity due to TiO and VO gases we term 
the "pM Class" planets, and those that are cooler, such that Ti and V are predominantly in solid condensates, 
we term "pL Class" planets. The optical spectra of pL Class planets are dominated by neutral atomic Na and K 
absorption. We calculate model atmospheres for these planets, including pressure-temperature profiles, spectra, 
and characteristic radiative time constants. Planets that have temperature inversions (hot stratospheres) of ^2000 
K and appear "anomalously" bright in the mid infrared at secondary eclipse, as was recently found for planets 
HD 149026b and HD 209458b, we term the pM Class. Molecular bands of TiO, VO, H 2 0, and CO will be seen 
in emission, rather than absorption. This class of planets absorbs incident flux and emits thermal flux from high 
in their atmospheres. Consequently, they will have large day/night temperature contrasts and negligible phase 
shifts between orbital phase and thermal emission light curves, because radiative timescales are much shorter 
than possible dynamical timescales. The pL Class planets absorb incident flux deeper in the atmosphere where 
atmospheric dynamics will more readily redistribute absorbed energy. This leads to cooler day sides, warmer 
' night sides, and larger phase shifts in thermal emission light curves. We briefly examine the transit radii for both 

classes of planets. The boundary between these classes is particularly dependent on the incident flux from the 
parent star, and less so on the temperature of the planet's internal adiabat (which depends on mass and age), and 
surface gravity. Around a Sun-like primary, for solar composition, this boundary likely occurs at ~0. 04-0. 05 AU, 
but uncertainties remain. We apply these results to pM Class transiting planets that are observable with the Spitzer 
Space Telescope, including HD 209458b, WASP- lb, TrES-3b, TrES-4b, HD 149026b, and others. The eccentric 
transiting planets HD 147506b and HD 17156b alternate between the classes during their orbits. Thermal emission 
in the optical from pM Class planets is significant red-ward of 400 nm, making these planets attractive targets for 
optical detection via Kepler, COROT, and from the ground. The difference in the observed day/night contrast 
between v Andromeda b (pM Class) and HD 189733b (pL Class) is naturally explained in this scenario. 
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1. INTRODUCTION a window into the compositi on of a planet and c ontains clues 

rm. , , . . . lt , . . „ ■ A t - , . to its formation history (e.g., Marlev et al. 2007). Of premier 

The blanket term hot Jupiter or even the additional term „, . , , — : — : , — . , 

, . T .. „, ,. ,■ .. c ,, . • . i • j ■ 4. i importance in this class of highly irradiated planets is how stel- 

very hot Jupiter belies the diversity of these highly irradiated , " , . „. , b / 

i 4. T3 i. i * i i i u v • 7 u • lar insolation affects the atmosphere, as this irradiation directly 

planets. Each planet likely has its own unique atmosphere, in- , , . r > j 

f . , , \- . ™ , 4. 4 c affects the atmospheric structure, temperatures, and chemistry, 

tenor structure, and accretion history. The relative amounts of , ,. r , ... , • , ■", 

c . A , .-, A ■ , . .,, a . the planet s cooling and contraction history, and even its stabil- 

refractory and volatile compounds in a planet will reflect the . F . b . ■" 

parent star abundances, nebula temperature, total disk mass, lo- ity against evaporation. 

.. flU , t , c ,. 4, ... , 4. c .. Since irradiation is perhaps the most important factor in de- 

cation of the planet s formation within the disk, duration of its ... , • • ,- , , 

c 4. , .4 , 4-4- \ T-u- termimng the atmospheric properties of these planets, we ex- 

rormation, and its subsequent migration (if any). This accre- . , b . ... ,,.,„„, ... 

, . . .,, . . . A . cc ■ 4 4 1 amine the insolation levels of the 23 known transiting planets, 

tion history will give rise to differences in core masses, total „. . , . & / „ 

, 1 4 if j u uj 4- We restrict ourselves to those planets more massive than Sat- 

heavy elements abundances, and atmospheric abundance ratios. ,, . , , F . , ,„ 

^- .,• ■ , .4 .. ■ 4, 4 c . , , urn, and hence for now exclude treatment of the hot Neptune 

Given this incredible complexity, it is worthwhile to first look _, T , . , . , , , , .. F , 

f , . , 4, . , . c , GJ 436b, which is by far the coolest known transiting planet, 

tor physical processes that may be common to groups of plan- „ m ' / ? 

Figure yj illustrates the stellar flux incident upon the planets as 

j , j j- c u 4 • a function of both planet mass (Figure ITfcf) and planet surface 

In addition to a mass and radius, one can further characterize . ,„ h-vn T , , , ■ 

, 4 u . , . .. 4 , ™ . , , • gravity (Figurellp). In these plots diamonds indicate transit- 

a planet by studying its atmosphere. The visible atmosphere is b J v b r 
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ing planets and triangles indicate other interesting hot Jupiters, 
for which Spitzer Space Telescope data exist, but which do not 
transit. 

The first known transiting planet, HD 209458b, is seen to 
be fairly representative of these planets in terms of incident 
flux. Planets OGLE-TR-56b and OGLE-TR-132b are some- 
what separate from the rest of the group because they receive 
the highest stellar irradiation. Both orbit their parent stars in 
less than 2 days and are prototype s of what has been called the 
class of "very hot Jupiters" ( Kon acki et al.ll2.Q03l : iBouchv et al.l 
2004) with orbital periods less than 3 days. However, orbital 
period is a poor discriminator between "very hot" and merely 
"hot," as HD 189733b clearly shows. Labeled a "very hot 
Jupiter" upon its d iscovery, due to its short 2.2 day period 
(Bou chv et alj|2005l) . HD 189733b actually receives a compar- 
atively modest amount of irradiation due to its relatively cool 
parent star. Therefore, perhaps a classification based on inci- 
dent flux, equilibrium temperature, or other attributes would 
be more appropriate. In this paper we argue that based on the 
examination of few physical processes that two classes of hot 
Jupiter atmospheres emerge with dramatically different spectra 
and day /night contrasts. Equilibrium chemistry, the depth to 
which incident flux will penetrate into a planet's atmosphere, 
and the radiative time constant as a function of pressure and 
temperature in the atmosphere all naturally define two classes 
these irradiated planets. 

Our work natur ally builds on the previous work of 
iHubenv et alj (2003) who first investigated the effects of TiO 
and VO opacity on close-in giant planet atmospheres as a func- 
tion of stellar irradiation. These authors computed optical and 
near infrared spectra of models with and without TiO/VO opac- 
ity. In general they found that models with TiO/VO opacity 
feature temperature inversions and molecular bands are seen 
in emissi on, rather than abso rption. Two key questions from 
the initial Hub env et alJd2003l) investigation were addressed but 
could not be definitely answered were: 1) if a relatively cold 
planetary interior would lead to Ti/V condensing out deep in 
the atmosphere regardless of incident flux, thereby removing 
gaseous TiO and VO, and, 2) if this condensation did not oc- 
cur, at what irradiation level would TiO/VO indeed be lost at 
the lower atmospheric temperatures found at smaller incident 
fluxes. 

Later iFortnev et all d2006bl) investigated model atmospheres 
of planet HD 149026b including TiO/VO opacity at various 
metallicities. Particular attention was paid to the temperature of 
the deep atmosphere P-T profiles (as derived from an evolution 
m odel) in relation to th e Ti/V condensation boundary. Similar 
to Huben vet al.l d2003l) . they found a temperature inversion due 
to absorption by TiO/VO and computed near and mid-infrared 
spectra that featured emissi on bands. Using the S pitzer In- 
fraRed Array Camera (IRAC) lHarrington et alJd2007h observed 
HD 149026b in secondary eclipse with Spitzer at 8 ;jin and de- 
rived a planet-to-star flux ratio consistent with a iFortnev et al.l 
(2006b) model with a temperature inversion d ue to TiO/VO 
opacity . At that point, looking at the wor k of Fortnev et al 



and especially IHubenv et alj d2003l) . [Harrington et al 



could have postulated that all objects more irradiated 
than HD 149026b may possess inversions due to TiO/VO opac- 
ity, but given the single-band detection of HD 149026b, cau- 
tion was in order. More recently, bas ed on the four-band de - 
tection of flux from H P 209458b by iKnutson et ail d2007al) . 
iBurrows et al. | d2007bl) find that a temperature inversion, po- 



tentialy due to TiO/VO opacity, is necessary to explain this 
planet's mid-infrared photometric data. Based on their new 
HD 209458b model and the previous modeling investigations 
these authors posit that planets warmer than HD 209458b may 
features inversions, while less irradiated objects such as HD 
189733b do not, and discuss that photochemical products and 
gaseous TiO/VO are potential absorbers which may lead to this 
dichotomy. 

We find, as has been previously shown, that those planets 
that are warmer than required for condensation of titanium (Ti) 
and vanadium (V)-bearing compounds will possess a temper- 
ature inversion at low pressure due to absorption of incident 
flux by TiO and VO, and will appear "anomalously" bright 
in secondary eclipse at mid-infrared wavelen gths. Thermal 
emission in the optical will be significant dHubenv et alj|2003l : 
lLopez-Morales & Seagerl2007l) . Furthermore, here we propose 
that these planets will have large day/night effective tempera- 
ture contrasts. We will term these very hot Jupiters the "pM 
Class," meaning gaseous TiO and VO are the prominent ab- 
sorbers of optical flux. The predictions of equilibrium chem- 
istry for these atmospheres are similar to dM stars, where ab- 
sorpti on by TiO, VO, H2O, and CO is prominent dLoddersI 
2002). Planets with temperatures below the condensation curve 
of Ti and V bearing compounds will have a gradually smaller 
mixing ratio of TiO a nd VO, leaving Na and K as the major op- 
tical opacity sources dBurrows. Marlev. & Sharpll2000h . along 
with H2O, and CO. We will term these planets the pL class, 
similar to the dL class of ultracool dwarfs. These planets will 
have relatively smaller secondary eclipse depths in the mid in- 
frared and significantly smaller day/night effective temperature 
contrasts. As discussed below, published Spitzer data are con- 
sistent with this picture. The boundary between these classes, 
at irradiation levels (and atmospheric temperatures) where Ti 
and V may be partially condensed is not yet well defined. 

In this paper we begin by discussing the observations to 
date. We then give an overview of our modeling methods and 
the predicted chemistry of Ti and V. We calculate pressure- 
temperature (P-T) profiles and spectra for models planets. For 
these model atmospheres we then analyze in detail the depo- 
sition of incident stellar flux and the emission of thermal flux, 
and go on to calculate characteristic radiative time constants 
for these atmospheres. We briefly examine transmission spec- 
tra before we apply our models to known highly irradiated giant 
planets. Before our discussion and conclusions we address is- 
sues of planetary classification. 

2. REVIEW OF SPITZER OBSERVATIONS 

2.1. Secondary Eclipses 

The Spitzer Space Telescope allows astronomers to mea- 
sure the thermal emission from ^3-30 /im from the highly 
irradiated atmospheres of these EGPs. This field has pro- 
gressed quickly from the first observations of secondary 
eclipses (when the plane t passes behind its p arent star) by 
ICharbonneau et all d2005l) for TrES-1 and by iDeming et al.l 
d2005l) for HD 209458b. Once it was clear what Spitzer's ca- 
pabilities were for these planets, additional observations came 
quickly. These incl uded secondar y eclipses for HP 1897 33b 
(Pe ming et al1l2006l) HP 1 49026b dHarrington et alj2007l) . GJ 
436b dPeming et aLlr2007l: iPemorv et alj|2007l) . and a wealth 
of new (yet to be published) data for ^10 other systems. In 
addition, mid infrared spectra were obtaine d for HP 209458b 
dRichardson et alj 120071: fSwain et al.ll2007l) and HP 189733b 



Two Classes of Hot Jupiter Atmospheres 



3 



dGrillmair et al]|2007l) . In particular, Ric hardson et al.l d2007l) 
claim evidence for e mission features in the mid-infrared spec- 
trum. Very recently Knuts on et al.l (|2007a) observed the sec- 
ondary eclipse of HD 209458b in all four InfraRed Array Cam- 
era (IRAC) bands. 

One can now begin to examine the secondary eclipse data 
for trends in planetary brightness temperatures. The brightness 
temperature, Tg is defined as the temperature necessary for a 
blackbody planet to emit the flux observed in the given ob - 
servational wavelength or band (Chamberl ain & Huntenf T987). 
For the published sec ondary eclipse data as of June 2007, 
lHarrington et al.l ((2007) show that the observed day side bright- 
ness temperatures for most transiting planets (TrES-1, HD 
189733b, HD 209458b) vary by ±20% from calculated plan- 
etary equilibrium temperatures (when making a uniform as- 
sumption for all planets regarding albedo and flux redistribu- 
tion). These differences are presumably due to wavelength de- 
pendent opacity and varying dynamical redistribution of the ab- 
sorbed flux. However, HD 149026b has a —45% higher bright- 
ness temperature than predicted from this relation, meaning that 
the temperature one measures at 8 /im is significantly higher 
than the planet's equilibrium temp erature. This lev el of infrared 
emission was predicted by iFortnev et alj d2006bl) for models 
of HD 149026b that included a hot stratosphere induced by 
absorption of stellar fl ux by TiO and VO. More recently, the 
iKnutson et all d2007iJ) observations of HD 209458b are broadly 
similar to the 1-band detection for HD 149026b. Relatively 
high brightness temperatures of 1500-1900 K were obse rved 
from 3.6 to 8 /im, compared to the Harring ton et al.l (|2007) pre- 
dicted equilibrium temperature of 1330 K. These observations 
have also been interpreted as b eing caused by a tem perature 
inversion (hot stratosphere) by iBurrows et al.l (l2007bl) . As we 
will discuss in §5.1, HD 209458b is in a temperature regime 
that is difficult to model, where Ti and V are expected to be 
partially condensed. 

2.2. Phase Curves 

The observations that best allow us to understand the atmo- 
spheric dynamics of these presumably tidally locked planets 
are those made as a function of orbital phase. Two approaches 
have been used to measure the thermal emission of these plan- 
ets away from secondary eclipse. The first is to periodically 
observe the infrared flux from the combined planet+star system 
at several times in the planet's orbit. The advantage of peri- 
odic short observations is that little telescope time is used. The 
disadvantage is that since the system is not monitored continu- 
ously, systematic uncertainties could be important and may be 
hard to correct for. There are published data for four planets 
that have been obse rved in this way: detec tions of phase vari- 
ation for v And b dHarrington et al.ll2006l) and HD 179949b 
dCowan et al.ll2007l) a nd upper limits for HD 209458b, 51 Peg 
b dCowan et alj|2007l) . Alternatively, one could monitor a sys- 
tem continuously over a significant fraction of a planet's orbital 
period, to eliminate any uncertainties induced b y having to re- 
visit the target. This method was employed by [Knutson et al.l 
d2007bl) to observe HD 189733b continuously for 33 hours, 
including the transit and secondary eclipse. Large day/night 
contrasts were found for pM Class planets v And b (which is 
consistent with a dark night side) and HD 179949b, whereas 
pL Class planet HD 189733b has a day /night 8 /im brightness 
temperature differe nce of only —240 K. We will return to the 
ICowan et al.l d2007l) upper limit for HD 209458b in §8. 



Two additional factors that affect these observations are the 
brightness of the parent star and the inclination of the plane- 
tary orbit. While the brightest planet-hosting stars allow for 
the largest flux measurements (e.g., v And b), the brightest hot 
Jupiter systems do not transit. A measurement of the secondary 
eclipse depth while obtaining a light curve is extremely valu- 
able, because the planet-to-star flux ratio at full planet illumi- 
nation is thus known, as is the planet's radius from transit obser- 
vations. The interpretation is then much more straightforward 
for the transiting systems. 

Among the current published data there appear to be two 
planets that are bright in secondary eclipse, HD 149026b and 
HD 209458b, compared to other transiting planets, and two 
planets that appears to have large day/night temperature differ- 
ences, v And b and HD 179949b. We will show that these plan- 
ets, HD 149026b, HD 209458b, v And b, and HD 179949b, are 
"pM Class" planets, while the less irradiated hot Jupiters (such 
as HD 189733b and TrES-1) are "pL Class." We will now turn 
to models of these classes of planets to examine how and why 
their atmospheres differ so strikingly. 

3. MODEL ATMOSPHERES 

3.1. Methods 

We have computed atmospheric pressure-temperature (P— 
T) profiles and spectra for several planets with a plane- 
parallel model atmosphere code that has been used for a 
variety of planetary and substellar objects. The code was 
first used to generate profiles and spectra for Titan's atmo- 
sphere by [McKay et al.l d 1989b . It was sig nificantly revised to 
model the atmospheres of brown dwarfs (Mar lev et al.l fl996: 
IBurrows etal] Il997t iMarley et al.l 120021) and irradiated giant 
planets llVIarlev & McKavl 19991 for Uranus). Recently it has 



been applied to L- and T-type brown dwarfs (Saumonetal 
| 2006[ l2007HCushing et al.l2007l) and ho t Jupiters dFortnev et a! 
|2005[ l2006bt IFortnev & Marlevll2007l) . It explicitly accounts 
for both incident radiation from the parent star and thermal ra- 
diation from the planet's atmosphere and interior. The radia- 
tive tr ansfer solution algorithm was developed by iToon et ail 
d 1989b . We model the impinging stellar flux from 0.26 to 6.0 
/im and the emitted thermal flux from 0.26 to 325 ixm. 

We use the elemental abundance data of Lodders (2003) 
and chemical equilibrium co mpositions are computed with th e 
CO NDOR code, followi ng lLodders & Feglevl d2002l 120061) . 
an d Lodders] dl99l|2002). We maintain a large and constantly 
updat ed opacity database that is described in iFreedman et al.l 
d2007l) . When including the opacity of clouds, such as 
Fe-metal and Mg-s i licates , we use the cloud model of 
lAckerman & Marlevl d2001l) . However, in our past work we 
have found onl y weak effects on P -T profiles and spectra due to 
cloud opacity dFortnev et al.| [2005). so we ignore cloud opacity 
here. However, the sequestering of elements into condensates, 
and their removal from the gas phase ("rainout") is always ac- 
counted for in the chemistry calculations. We note that day 
sides of the strongly irradiated pM Class planets are too warm 
for Fe-metal and Mg-silicates condensates to form. 

3.2. Very Hot Jupiters and TiO/VO Chemistry 

It is clear that the abund ances to TiO and VO gases is impor- 
tant in these atmospheres. Hube nv et al.l ((2003) first discussed 
how understanding the "cold trap" phenomenon may be signif- 
icant in understanding these abundances. If a given P-T profile 
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crosses a condensation curve in two corresponding altitude lev- 
els, the condensed species is expected to eventually mix down 
to the highest pressure condensation point, where the cloud re- 
mains confined due to the planet's gravitational field. This pro- 
cess is responsible for the extremely low water abundance in 
the Earth's stratosphere. It can also be seen in the atmospheres 
Jupiter and Saturn, where the ammonia ice clouds are confined 
to a pressure of several bars, although both these planets exhibit 
stratospheres, such that their warm upper atmospheres pass the 
ammonia condensation curve again at millibar pressures. For 
the highly irradiated planets, the relevant condensates are those 
that remove gaseous TiO and VO, and sequester Ti and V into 
solid condensates at pressures of tens to hundreds of bars, far 
below the visible atmosphere. 

The cold trap phenomenon constitutes a departure from 
chemical equilibrium that cannot be easily accounted for in 
the pre-tabula t ed ch e mical equilibri u m abun d ances used by 
Hub enyet alJ (120031) . IFortnev et al] (l2006bl) . iBurrows et all. 
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d2007bl) . and here as well. Our chemical abundances and opac- 



ities are pre-tabulated in P-T space and the atmosphere code 
interpolates in these abundances as it converges to a solution. 
The abundances determined for any one pressure level of the P- 
T profile are not cognizant of abundances of other levels of the 
profile, although condensation and settling of species is always 
properly accounted for. In this case a tabulated TiO abundance 
at a given P-T point at which, in equilibrium, TiO would be in 
the gaseous state (warmer than required for Ti condensation), 
may not be correct. If the atmospheric P-T profile intersects 
the condensation curve, the atmosphere becomes depleted in 
TiO above the cloud. We do not treat the cold trap here. In 
practice, we use two different opacity databases: one with TiO 
and VO removed at P<10 bars, which simulates the removal 
of Ti and V into clouds, and one in which gases TiO and VO 
remain as calculated by equilibrium chemistry (which does not 
include depletion from a cold trap). 

A full discussion of titanium and vanadium chemistry in 
the cont e xt of M- and L-dwarf atmospheres can be found in 
lLoddersI d2002l) . Much of that discussion pertains to the atmo- 
spheres of highly irradiated plan ets as well. Th e chemistry is 
complex. For instance ILodders & Feglevl (12006) (using the up- 
dated Lodders, 2003 abundances) find the first Ti condensate 
will not necessarily be CaTi03. For solar metallicity the first 
condensate is TiN if P > 30 bars, Ca 3 Ti 2 7 if 5 < P < 30 bars, 
Ca4Ti 3 OH) if 0.03 < P < 5 bars, and CaTi0 3 if P < 0.03 bars. 
These four condensates are the initial condensates as a function 
of total pressure and their condensation temperatures define the 
Ti-condensatio n curve in Fig ure[2] Another important point is 
that, following lLoddersI d2002l) we assume that vanadium con- 
denses into solid s olution with Ti-bearing c o ndensates, as is 
found in meteorites dKornacki & Feglevll986l) . lLoddersI (|2002) 
find this condensation sequence is also co nsistent with ob- 
serve d spectra at the M- to L-dwarf transition ( Kirkpa trick et alj 
119991) . In the absence of this effect V would not condense until 
^200 K cooler temperatures are rea ched and solid VO forms , 
as in the ch e mistry calc ulations of Burrows & Shard ([1999), 
lAllard et all d2001l) . and ISharp & Burrows! d2007l) Choices 
made in the calculation of the Ti condensation curve (e.g., num- 
ber of potential Ti-bearing compounds included in the calcula- 
tions) and the V condensation curve (e.g. whether V enters into 
Ti-bearing condensates or condenses as solid VO) will cause 
some shift in the position of these condensation curves, impor- 
tant boundaries between the pL and pM Classes. 



3.3. Calculating Pressure-Temperature Profiles 

When modeling the atmospheres of irradiated objects with 
a one-dimensional, plane parallel atmosphere code it is neces- 
sary to weigh the stellar flux by a geometric factor if one is 
computing a profile for day-side average or planet-wide aver- 
age conditions. Decriptions of this issue in solar system at- 
mos phere modeling can b e found in I Appleby & Hoganl (Q984) 
and iMcKav et alj d!989l). In the context of highly irradi- 
ated EGPs. IBurrows et all (120041 120061 l2007al). iBarman et al.1 
d2005h . lSeager et alJd2005l) , and lFortnev & Marlevld2007l) have 
all discussed this issue in some detail. All are approxima- 
tions for atmospheric structures that are surely complex for 
these highly irradiated atmospheres. In this paper, we make 
a straightforward choice to multiply the incident flux by a fac- 
tor, /=0.5 for a day-side average (since the atmosphere radiates 
over 2ir steradians but intercepts flux over only ir radians), and 
assume normal incidence of flux, /i=l. T his differs slightly 
from choices in our previous papers (see IFortnev & Marlevl 
2007). All profiles shown in this paper are relevant for the irra- 
diated planetary day side. 

In Figure|2]we show P-T profiles as a function of distance 
from the Sun. At 0.025 and 0.03 AU the profiles are every- 
where warmer that the condensation curve for a solar abun- 
dance of Ti and V. Essentially all Ti and V is expected to be in 
TiO and VO. By 0.035 AU, where the incident flux is weaker, 
the profile crosses this condensation curve, labeled "IX Ti/V- 
Cond." To the left of this curve, the TiO/VO abundances fall 
exponentially with decreasing temperature. To the right, the 
TiO/VO abundances are essentially constant. This particular P- 
T profile is oddly shaped since it samples a region of P-T space 
where the opacity drops rapidly with temperature. It crosses 
back to the right of the condensation curve, where it again finds 
the "full" abundance of TiO/VO. This is an example of a pro- 
file that violates the cold-trap (which we do not model), as at 
lower pressures TiO/VO abundances should be reduced, given 
the condensation of Ti and V below. The dotted curve labeled 
"90% Ti-Cond" shows where 90% of Ti has been lost to a solid 
condensate (Lee & Lodders, in prep). Figure[3]is a plot of opti- 
cal and near infrared opacity at 1 mbar and shows that TiO/VO 
are still major optical opacity sources even after condensation 
has begun. At 1700 K, before Ti/V condensation, opacity (here 
log(K) in cmV 1 ) in V-band is —0.5. By 1600 K, 33 K cooler 
than the "90% Ti-Cond" curve, this opacity is still an order of 
magnitude larger than the opacity at 1400 K, by which point the 
strong TiO/VO bands have given way to alkali lines and water 
bands. 

That TiO/VO gradually wane in importance is c onsistent 
with observations of dwarfs at the M-to-L transition (Lodders 
2002). The mixing ratios of these species do not drop to zero 
at the solar abundances condensation curve. For model atmo- 
spheres more distant from the parent star (0.05 AU and beyond), 
the abundance of TiO/VO continues to fall exponentially with 
temperature, such that in the upper atmosphere it is no longer a 
major opacity source. 

We can examine how the derived P-T profiles for the pM 
Class planets vary as a function of 7i nt , a planet's intrinsic ef- 
fective temperature in the absence of irradiation, and as a func- 
tion of surface gravity. This is shown in Figure [4] Lower grav- 
ity planets have lower pressure photospheres, and consequently 
have warmer interiors, all else being equal, which has long been 
known from basic theory. As can be seen in Figure[TJ?, the low- 
est gravity planets cluster around g=9 m s" 2 , while the high- 
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est gravity planets, with the exception of HD 147506b, cluster 
around g=25 m s" 2 , similar to Jupiter. Both of these values are 
a factor of 2/3 removed from our baseline g=15 m s~ 2 . The ef- 
fect of gravity on the P-T profiles for the pM Class planets are 
relatively modest. However, the effect of 7i nt is potentially im- 
portant. We show profiles at 7i nt values of 150, 200, and 250 K. 
Generally for adiabatic interiors, the higher the 7j nt , the warmer 
the interior, and the larger the radius of the planet. Independent 
of 7i n t, which could in principle be affected by barr iers to con- 
vection (IStevensonll985HChabrier & Baraffe 2007), the higher 
the entropy of the interior adiabat, the warmer the atmospheric 
adiabat, the larger the radius of the planet, and the less likely 
the deep atmosphere profile at 10 to 1000 bar will cross the 
Ti/V condensation boundary. Over Gyr-timescales it may be 
possible for a planet's interior to cool enough such that signifi- 
cant Ti/V condensation could occur at depth, and a planet could 
transition from pM to pL Class. 

In Figure Q] the gray shaded area is meant to illustrate a tran- 
sition region between these classes of planets, where TiO and 
VO are waning in importance. Given the lack of data for the 
systems in this proposed transition region, the exact location 
and height of this region is merely suggestive. Analysis of the 
current secondary eclipse data for HD 189733b dDem ing et al. 
120061: iGrillmair et alj2007HKnutson et alj2007b|) fits well with 
a pL Cl ass model (fFortney & Marlevll2007l) . and Knutson et al.l 
d2007bl) directly measured a small day-night contrast, suggest- 
ing that TiO/VO opacity has waned considerably at this incident 
flux level and there is no temperature inversion. We next turn 
to the upper boundary. 

At the high incident fluxes encounter ed by HD 149026b , 
the temperature inversion predicted by Fortne vet alj (|2()06b) 
neatly explains the lar ge planet-to-star flux ratio at 8 /xm 
dHarrington et al.l [20071) . At lower incident fluxes, the large 
day-night temperature contrasts implied for v And b and HD 
179949 imply temperature inversions as we'll show in §4.2 (and 
pM classificatio n) at this flux level. Intr iguingly, at a still lower 
irradiation level. iBurrows et all d2007bl) find that a strong tem- 
perature inversion is needed for HD 209458b, which we also 
address in §5.1, which indicates that inversions extend down to 
at least this incident fluxs level. This inversion persists even 
though it is clear from Figure[2]that Ti/V condensation has be- 
gun in our ID model at this irradiation level (although is has 
not for the more irra diated v And b , as we show in §5.2). Simi- 
lar to the findings of lBurrows et al. (2007b), we can only obtain 
a temperature inversion for HD 209458b induced by TiO/VO 
opacity if their abundances are not diminished by the simple 
cold trap description discussed above. In §5.1 and §7.1 we dis- 
cuss a few avenues to towards understanding if TiO/V O remains 
the prominent stratospheric absorber down to this incident flux 
level. 

The transition area as plotted in Figure[T]is thus given from 
mix of the observations and theory available to date. Addi- 
tional complications will be effects of atmospheric metallicity 
on Ti/V condensation, the treatment of the incident stellar flux, 
planet surface gravity, and the temperature of the interior adi- 
abat. For instance, it is certainly possible that a planet could 
appear marginally on to the hot side of this divide, but with a 
very cool interior (due to low mass and/or old age) would have 
Ti/V condensed such that it is a pL Class member. With better 
knowledge of the incident flux deposition and Ti/V condensa- 
tion, it may eventually be possible to "take the temperature of 
the interior adiabat" at pressures of 10 to 100 bars, which could 



shed light on the interior structure of some of these planets. 

4. RESULTS AND PREDICTIONS 

4. 1 . Absorption and Emission of Flux 

It is worthwhile to examine in some detail how the atmo- 
spheres of the pM Class and pL Class planets differ in the ab- 
sorption of stellar flux as a function of atmospheric pressure. 
One can plot the layer net fluxes for several pressures as a func- 
tion of wavelength to illustrate the main atomic and molecu- 
lar absorbers in the atmosphere. Since these models are in ra- 
diative equilibrium, the integrated absorbed flux within a layer 
should equal the integrated emitted flux from that same layer. 
We can plot the net and integrated thermal flux to examine at 
what wavelengths, and which species, lead to atmospheric heat- 
ing and cooling. A series of plots are shown in Figures [5] and[6] 
for a pL Class planet at 0.05 AU, and Figures[7]and[8] for a pM 
Class planet at 0.03 AU. The left ordinate of these plots shows 
the quantity erg s" 1 g" 1 fjanT 1 , which is the power per wave- 
length interval, per gram of atmosphere, absorbed or emitted in 
a slab of atmosphere, though we'll maintain the term "flux" for 
convenience. 

In Figure[5]five panels are shown, from 0.45 mbar to 4.1 bar, 
each separated by approximately an order of magnitude in pres- 
sure. Prominent in the top two panels are absorption of incident 
flux by neutral atomic alkalis Na (at 0.59 /jm) and K (at 0.77 
fim) as well as strong H2O bands in the near infrared. While Na 
and K lead to strong absorption in optical wav elengths, caus- 
ing extremely dark atmo spheres in the optical dSudarskv et al.l 
2000l lSeager et al.ll2000l) . at 0.45 and 4.4 mbar, respectively, 4 
and 1 .5 times more incident stellar is actually absorbed by H2O 
than by alkalis. The third panel, 43 mbar, shows that essen- 
tially all incident flux at the wavelengths of the Na and K line 
cores have been absorbed, and the alkali line wings are now as 
important as the H2O bands in absorption. By 420 mbar, the 
wavelengths corresponding to the broad Na and K lines have 
already become optically thick, leaving little flux left. About 
~10 3 less flux is absorbed at 4.1 bar than 0.45 mbar. The pres- 
sure is now in the nearly isothermal region of the P-T profile, as 
shown in Figure [2] In summary, Na, K, and H2O are the major 
absorbers and heating in the millibar region of the atmosphere 
amounts to ~ 10 6 erg g -1 s -1 . 

We can examine Figure |6]to understand the cooling of the at- 
mosphere as these same pressure levels. This figure has a much 
larger wavelength range that is shown on a log scale. In the top 
three panels, from 0.45 to 43 mbar, it is predominantly the near 
and mid infrared water bands that allow the atmosphere to cool 
to space. Between ^0.8-1.8 /im, positive layer net fluxes occur 
due to absorption of flux from above, by H2O. Deeper in the 
atmosphere, at higher temperatures, the local Planck function 
at a given layer overlaps shorter wavelength water bands from 
~l-2 /im that are primarily responsible for cooling. 

In Figure [7] we can clearly see how this pM Class planet dif- 
fers from the pL Class. The incident stellar flux is higher and 
more of this flux is deposited into the 0.45 mbar layer of this 
model, as compared to that shown in Figure[5] This deposi- 
tion occurs almost entirely through absorption by TiO and VO 
bands in the optical and near infrared. Water vapor still absorbs 
flux, but its effect is swamped by the TiO/VO absorption. This 
continues at 4.4 mbar, but already by 43 mbar, the vast ma- 
jority of the incident flux has been absorbed, especially in the 
optical wavelengths. By 420 mbar there is little flux left in the 
optical and by 4. 1 bar there is little flux left at any wavelength; 
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again this pressure is in the nearly isothermal region shown in 
Figure |2 

The top panel of Figure[8] shows that, since the atmosphere 
has absorbed a significant amount of flux at low pressures, com- 
pared to the pL Class model, that the atmosphere must reach a 
much higher temperature to be able to adequately radiate away 
this energy. The local Planck function moves significantly blue- 
ward, so that it can radiate by means of the strong H2O near in- 
frared bands as well as the optical bands of TiO and VO. At 4.4 
and 43 mbar, the atmosphere is cooler, and radiation by H2O 
bands is more important. At 43 mbar some heating from above 
does occur via TiO/VO absorption. At 420 mbar and 4. 1 bar, 
the local temperature has again increased, leading to radiation 
again most prominently by the near infrared H2O bands. 

4.2. Radiative Time Constants 

It is clear that a key difference between the atmospheres of 
the pL Class planets and pM Class planets is the pressures at 
which the absorption and emission of flux occurs. This can 
be shown in a more straightforward manner by means of the 
brightness temperature, T%. One can then examine the P-T pro- 
file to find the pressure that corresponds to a given Ts, which 
is the characteristic atmospheric pressure for this thermal emis- 
sion. This quantity is plotted as a function of wavelength in 
Figure [9] In general thermal emission arises from a pressure 
level roughly an order of magnitude greater in a pL Class atmo- 
sphere than in a pM Class atmosphere, due to the higher opac- 
ity in the hotter pM atmospheres, both at optical wavelengths, 
due to TiO/VO, and infrared wavelengths, due to increased H2O 
opacity and H2 collision induced absorption. 

What effect this may have on the dynamical redistribution of 
energy in a planetary atmosphere can be considered after cal- 
culation of the radiative time constant, r, a( j. In the Newtonian 
cooling approximation a temperature disturbance relaxes expo- 
nentially toward ra diative equilibrium with a characte ristic time 
constant r rad (e.g.. lGoodv & Yung|[l9l^lSalbvl [i996). Atpho- 
tospheric pressures this value can be approximated by 
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where a is the Ste fan-Boltzmann constant an d cp is the spe- 
cific heat capacity dShow man & Guilloti r2002). However, this 
quantity can be derived anywhere in the radiative atmosphere 
via the following formulation. We first obtain a P-T profile so- 
lution that is in radiative-convective equilibrium. We can then 
take this profile, and in a given layer include a small (^10 K) 
thermal perturbation, AT. With this AT in place we perform 
a radiative transfer calculation, and additionally calculate the 
flux divergence dF/dz, which is zero in radiative equilibrium 
(where z is the height). For a function /(f) a time constant will 
be given by / x (df/dt)~ l . We can then solve for r ra d as 

where p is the mass density. The quantity (dF jdz) /{pep) is the 
heating/cooling rate in K s . By varying the location of the 
AT perturbation with height, one can calculate r ra d as a func- 
tion of pressure in the atmosphere, for a given P-T profile. In 
practice we find values o f r ra( j that are very similar to those de- 
rived bv llroet all d2005l) for HD 209458b. Interestingly, these 
authors performed a much different calculation that involved 
the input of Gaussian temperature fluctuations and the direct 
calculation of the time necessary for their code to relax back to 



radiative equilibrium. A larger exploration of r ra d over a wider 
range of phase space, with additional discussion, will be found 
in A. P. Showman et al. (in prep.). 

In Figure[10] we show our calculated r ra d as a function of 
pressure alongside the radiative-convective equilibrium profiles 
from which they were generated. Although these time constants 
are nearly equivalent in the dense lower atmosphere, they are 
significantly different in the thinner upper atmosphere that one 
is sensitive to from mid-infrared observations. 

The right y-axis in Figure [9] is cast in terms of the radia- 
tive time constant appropriate for a given pressure in the at- 
mosphere. Note that this axis is not linear, but each major tick 
mark is labeled with the r ra d that is appropriate for the major 
tick mark at the left. While those for the pM Class range from 
only 10 3 — 10 4 s, those for the pL Class range from 1 4 — 1 5 
seconds. This ~ 10 x difference in timescale is a consequence 
of both the hotter temperatures and lower pressures of the pM 
Class photospheres, as can be understood from EquationQ] 

One can also define an advective time scale, a characteristic 
time for moving atmospheric gas a given planetary distance. A 
common definition is 

7"advec — "7T 5 (3) 

where U is the wind sp e ed and R p i s the planet radius 
(I Showman & Guillotl 120021 ISeager et al.1 l2005h . If one sets 
Tadvec = 7"rad an d sets = l-Rjs one can derive the wind speed U 
that would be necessary to advect atmospheric gas before a r ra d 
has elapsed. This is also shown in Figure|9] via the gray bars on 
the right side. Since depth to which one "sees" is wavelength 
dependent, the apparent efficiency of atmospheric dynamics in 
redistributing energy in the planet's atmosphere will depend on 
the wavelength of observations. For wavelengths where the at- 
mospheric opacity is low, one sees flux emerging from rela- 
tively deep layers. The planet will appear more homogenized 
( Seager et al.l l2005). as was sh own explicitly by iFqrtnev et al.l 



(2006a) for light curves for the lCooper & Showm an (20061) 3D 
dynamical model. Therefore the shape and times of maxima 
and minima in thermal emission light curves will be a function 
of wavelength. 

Examination of Figure|9] shows that for a pM Class planet 
winds speeds of dozens of km s -1 would be necessary for ad- 
vection to dominate over cooling/heating ! For a pL Class planet 
more modest winds speeds of several km s" 1 are needed. Sev- 
eral km s" 1 winds, which are similar to the sound speed, are 
in line with predictions of a number of 2D and 3D dynam- 
ical models for hot Jupiter atmospheres dShowman & Guillotl 
2002 1 ICooper & Showmanl 120051; iLangton & Laughlinl I2007F 
Dobb s-Dixon & Linl [2007h . These calculated wind speeds 
are far below those needed to advect air before it radiatively 
heats/cools in a pM Class atmosphere. 

Winds will not be able to redistribute energy at the photo- 
spheres of pM Class planets. The atmospheres of pM Class 
planets likely appear as one would expect from radiative equi- 
librium considerations: the hottest part of the atmosphere is at 
the substellar point, and the atmosphere is cooler as one moves 
toward the planet's limb. Deviations from this prediction may 
occur from 1-2 /im where opacity windows allow for obser- 
vation of flux from relatively higher pressures. The night-side 
temperature will be relatively cold and will be set by the effi- 
ciency of energy redistribution at depth, as well as the intrinsic 
flux from the interior of the planet. In the absence of a large 
intrinsic flux driven by an interior energy source, such as tides, 
redistribution may well swamp the intrinsic flux. For a pL Class 
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atmosphere, there will be a much more complex interplay be- 
tween radiation and dynamics; energy redistribution will more 
readily lead to a planetary hot spot being blown downwind. The 
location of the hot spot will itself be wavelength dependent, 
and there will increased energy t ransfer to the night side at 
photospheric pressures. Recently iDobbs-Dixon & Lin] ((2007) 
have shown that larger day-night temperature differences are 
expected with increased atmospheric opacity, but they models 
lack motivation for their opacity choices. Dynamical models 
that do not include a realistic treatment of opacities and radia- 
tive transfer will miss the important differences between these 
two classes of planets. 

4.3. Characterization Via Transit Observations 

Primary transit observations for pL Class planets HD 
209458b and HP 189733b have yielded detec tions of sodium 
(Charbonneau et al. 2002) and water vapor dBarmanl [2007; 
iTinetti et al .1120071) . Since water vapor is present in abundance 
in pM Class planets, it should be detectable as well. Observa- 
tions in the optical would show strong absorption due to TiO 
and VO across a bro ad wavelength range , which will lead to 
observed transit rad ii dHubbard et : alJl200ltlBaraffe et al.|[2003t 
Burr ows et al.l2003l) that correspond to even lower atmospheric 
pressure (higher in the atmosphere) that those forpL Class plan- 
ets. 

We perform a simple calculation to illustrate this point. 
From Figure|9] the approximate pressure levels for thermal 
emission are known. After calculating the distances between 
the wavelength-dependent photospheric pressures to a refer- 
en ce pressure of 1 bar, we can follow the implementation 
of iBurrows et alj d2007al) to derive the (even larger) radius at 
which the slant optical depth reaches unity. This is the radius 
one measures (at a give n wavelength or waveb and) during a 
transit observation. From lBurrows et aTl (1200 7a). the extension 
dR from the photospheric radius to the transit radius is approx- 
imately 

dR = H\n\j-^, (4) 

where H is the characteristic atmospheric scale height and R p is 
the planet radius. We plot the approximate transit radius for the 
model pL Class planet and pM Class planet in Figure [TT] Here 
the 1 bar radius is 1.20 Rj and the characteristic T for the scale 
height calculation for the pL and pM Class planets are 1070 
K and 2000 K, respectively. Since our derived P-T profiles do 
not reach to microbar pressures, a more detailed transit radius 
calculation utilizing the precise run of temperature and opac- 
ity vs. pressure (see [Hubbard et aTl l200U IFortnev et ail 12003) 
will be investigated in future work. Nonetheless, the important 
points are clear. Due predominantly to the hotter atmosphere of 
pM Class planets, their at mospheric extension is larger. Ad- 
ditionally, as suggested by IBurrows et a D (|2007b) the strong 
opacity of TiO and VO leads to a larger radius in a wide optical 
band, compared pL Class planets. Figure [TT] shows a differ- 
ence of ^5% at optical wavelengths. This difference in radius 
could be even further enhanced by temperatures in excess of 
the 2000 K used in the calculation. This plot should only be 
considered suggestive, as a more detailed calculation would in- 
volve different chemical abundances on the day and night side 
of the planetary limb. Cooler limb temperatures may not allow 
for gaseous TiO and VO. 

Detailed secondary eclipse observations may shed additional 
light on the structure of these atmospheres. As was shown in 



IFortnev et all d2006bl) for HD 149026b, planets with tempera- 
ture inversions will show limb brightening, rather than darken- 
ing, whic h may eventually be detectable via secondar y eclipse 
mapping dWilliams et al.1 12006; Rau scher et alj|2007bl) . 

5. APPLICATION TO KNOWN PLANETS 
5.1. HD 209458b at Secondary Eclipse 

Recently Knut son et alj d2007al) published observations of 
the secondary eclipse of planet HD 209458b. These observa- 
tions utilized Spitzer IRAC, and constitute the first published 
observations for a transiting planet across all four IRAC bands. 
These observed planet-to-star flux ratios were interpreted as 
being caused by a temper ature inversion (hot stratosphere) on 
the day side of the planet dBurrows et alj|2007bt iKnutson et al.1 
2007a). The reasons for this interpretation include a large 
brightness temperatures relative to its expected day-side T e s, 
especially at 4.5 and 5.8 /j,m, and a large flux ratio at 4.5 
/mi (relative to 3.6 /im), which is enhan ced due to strong CO 
and H2O emission fe atures from 4-6 /im dFortnev et alj |2006b; 
IBurrows et aT]|2007bl) . In our terminology, HD 209458b joined 
HD 149026b as a pM class planet. 

In FigureQjh we plot P-T profiles for HD 209458b, with 
and wi thout the opacity of TiO and VO. As was previously 
shown dHubenv et alj2003l:rFortnev et alj2006bllBurrows et al.l 
2007bl) absorption of stellar flux by TiO/VO leads to a hot- 
ter upper atmosphere and cooler lower atmosphere, for a given 
planet. The profile that includes TiO/VO opacity (red) is cooler 
than the condensation curve for a solar-abundance mixture of 
Ti/V. For this profile, given our equilibrium chemistry table, 
the abundance of TiO gas drops even below 10% of a solar Ti 
abundance, but then increases back to solar in the hot upper 
atmosphere. A proper treatment of the cold trap could leave 
only a fraction of a solar abundance of Ti/V in TiO/VO in the 
upper atmosphere. That being said, a hot stratosphere similar 
to one that we derive is nessary to explain the planet's spec- 
trum. We think that opacity due to TiO/VO is a natural ex- 
planation, but more work needs to be done to understand the 
abundances of these molecules at low pressure in these dy- 
namic atmospheres , such as a non-equilib rium cloud conden- 
sation model (see lHelling & Woitk e 2006, and §7.1). In green 
is a similar profile that utilizes a different me thod of treating the 
angul ar dependence of the incident flux (see lFortnev & M arlev 
120071) . In blue is the day-side profile computed with TiO/VO 
opacity removed. 

Figure \l2b shows the spectra for these three models. The 
flux levels at 4.5 and 5.8 /im are so large that they can only be 
caused by a temperature inversion. This is similar to wha t was 
found at 8 /im for HD 149026b bv lHarrington et alj d2007l) . In 
the red model the temperature gradient at several millibars is 
not steep, so the emission features are quite weak. Emission 
features are more prominent in the green model. Both models 
with a temperature inver sion also fit the 2.2 /im constraint from 
iRichardson et alj d2003l) 

IBurrows et al.l d2007b) present a nice fit to the IKnutson et al.l 
d2007al) observations. With their models they explore two ad- 
ditional parameters not found in our models. Most importantly, 
since they include a sink of energy on the day side, this cools 
their HD 209458b atmosphere significantly (down to 1000 K) at 
the pressures responsible for emission in the 3.6/im band, lead- 
ing to much less flux in this band. This enables them to simul- 
taneously fit the 3.6 and 4.5 /im band observations. In addition, 
the opacity source that causes their temperature inversion has an 
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optical depth for absorption that is similar to what they predict 
for Ti O/VO, but it is not TiO /VO specifically. Similar to our 
model, Burro ws et alj d2007bl) have trouble simultaneously fit- 
ting the 5.8 /im and 8.0 /im observations. Due to strong opacity 
due to water, emission at 24 /mi emerges from high in the at- 
mosphere (see Figure[9]). In our models , this pressure is within 
the 2000+ K hot stratosphere, whereas iBurrows et all d2007b) 
finds a cooler upper atmosphere (~1500 K) above their inver- 
sion. In a future work we will explore P-T profiles that ex- 
tend to lower atmospheric pressure, to better understand how 
this may effect our derived profiles. The most striking aspect 
of the lKnutson et al.l ((2007a) observations is the relatively cool 
temperature needed to fit the 3.6 and 8.0 /im bands, relative to 
those at 4.5 and 5.8 /im. Although our day side fits shows room 
for improvement, it is exciting that this planet falls into a class 
with HD 149026b. The details of absorption of incident flux by 
atmospheres partially depleted in TiO/VO due to condensation 
must be worked out in more detail before it is clear to what de- 
gree the TiO/VO abundances in the atmosphere of HD 209458b 
are anomalously high. 

HD 209458b is a long-term ta rget for the detect ion of opti- 
cal flux from the MOST satellite. iRowe et a l. (2007) find a new 
lcr upper limit to the planet-to-star flux ratio. Figure[T2"t> shows 
that both models of HD 209458b, with and without a tempera- 
ture inversion, are consistent with this limit. For the pL Class 
model, this is due to strong optical pressure-broaden ed lines of 
Na and K, as has been ex pected for some time ( Sud arskv et al.l 
l2000t ISeager et alJl2000h . For the pM Class model TiO/VO 
opacity leads to less scattering of stellar flux, and TiO/VO band 
emission is not bright enough to lead to a detection. It is a 
bit disappointing that MOST bandpass (^400-680 nm) happens 
to coincide with a minimum in planetary flux. However, the 
COROT red bandpass and Kepler band extend to redder wave- 
lengths (~ 1000 nm and ~ 850 nm respectively) such that the 
detection of planetary optical flux by these satellites for pM and 
pL Class planets should be achievable. 

The opacity of TiO/VO leads to yet another explanation for 
the weaker than expected sodium absorption feature seen by 
ICharbonneau et af] (12002b. Although condensates are still an 
attractive option dFortney et al.l2003UFortneyl20'05l) . if TiO/VO 
are found on the day side of the planet's limb, opacity due 
to these molecules from ~ 400 to 1000 nm could swamp the 
expected Na absorption feature at 589 nm. If the opacity on 
the limb is everywhere larger in the optical this would lead 
to weaker, narrower Na absorption peak than expected by the 
first-generation mod els dSeager & Sasselovl2000l:lBrownl200lt 
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This effect can be seen in calculations of 
d2002l) . whose description of atmospheric chem- 
istry at that time left addition absorbers of optical flux in the 
upper atmosphere. It is not clear if TiO/VO would be found in 
abundance at the planet's limb, but this is an avenue worthy of 
further exploration. 

5.2. Other Planets 

We can examine the atmospheres of likely pM Class plan- 
ets from Figure Q] Since the OGLE planets will be difficult 
to detect in the infrared, we show P-T profiles for W ASP- lb 
dCollier Cameron et al.ll2007t ICharbonn eau et al 1l2007l). TrES - 
4b dMandushev et al.ll2007l). TrES-3b dO' Donovan et alJl2007h . 
and v And b dButler et all 119971: iFuhrmann et al.l 1 19981) in 
Figure Qj] For these profiles the temperature of the adiabatic 
deep atmosphere at tens to hundred of bars (the top of the inte- 



rior adiabat) has been chosen to be consistent with the measured 
planetary radius, assuming a J upiter/Saturn like interior abun- 
dance of ice/rock (25 M®, see[Saumon & Guillotll2004l) . The 
chosen adiabats do not rely on any specific evolution model — 
for a given planet mass and composition only one particular 
inter nal adiabat will mat ch the planet's measured radius (see, 
e.g., iHubbard et al . 2001, Figure 1). For the large-radii transit- 
ing planets, deep radiative zones extending down to ~1 kbar 
are not consistent with the relatively warm interiors implied by 
these radii. For v And b, which has a mass similar to that of 
HD 209458b, radius and surface gravity values for HD 209458b 
were used. The temperature of the upper atmosphere of these 
models are a function of the incident stellar flux, which de- 
creases in magnitude from WASP- lb to v And b. We note 
that for these hotter pM Class planets, the day-side profiles are 
everywhere warmer than that required for Ti/V condensation, 
such that opacity due to TiO/VO should not be considered at 
all anomalous. For comparison we also pl ot a pL Class model 
for potential tran sition object HAT-P-lb dBakos et al.ll2007bt 
IWinnetalJl2007h . which transits one member of a main se- 
quence G star binary system. We highlight it here because 
the binary nature of the system makes it a good candidate for 
follow-up observations that rely on relative photometry. 

We can examine the computed planet-to-star flux ratios for 
these models. These are shown in Figure[l4] using the same 
color scheme. Most striking is TrES-3b, whose large planet- 
to-star radius ratio leads to large flux ratios reaching 0.01 at 
longer mid-IR wavelengths, as seen in Figure [14b. The pM 
Class models have fairly shallow temperature gradients in the 
hot stratosphere, leading to relatively week emission features in 
the spectra. A comparison with HAT-P-lb shows that weak wa- 
ter emission features are seen in the pM Class planets at wave- 
lengths where HAT-P-lb shows water absorption bands. The 
feature at 4.5 /mi is due predominantly to absorption in the 
stellar atmosphere, along with planetary emission due to CO. 
The incident stellar spectra are tailored for ea ch planet's parent 
star an d orbital distance, and are taken from Haus childt et all 
(fl999l 

Figure \14b shows ratios in the optical and near in- 
frared. As has been pointed out previously, for these 
hot atmospheres the rmal emission con tr ibutes significantly 
to the optical flux dFortnevetal.1 120051; ISeager et all [2005; 
lLopez-Morales & Seagerll2007l) . These pM Class planets will 
be excellent targets for red-optical secondary eclipse observa- 
tions, as the planet-to-star flux ratio reaches ~ 10~ 4 . It is im- 
portant to note that this optical flux is thermal emission, not 
scattered incident flux. For the TrES-4b model {red), the ver- 
tical bar indicates that by 0.45 /im, thermal flux is already 
10 times greater than reflected light. Thermal flux dominates 
at all longer optical wavelengths. Even for pL Class model 
for HAT-P-lb, thermal emission is 10 times greater than re- 
flected light at 0.73 /im. If one wanted to observe true "re- 
flected light" from these atmospheres, the blue spectral regions 
should be targeted. In a future work we will examine geo- 
metric albedos and lightcurves for these planets to understand 
current and future data from spacecraft like MOST, Kepler, 
and COflOr. lLopez-Morales & Seagerl d2007l) have recently in- 
vestigated prospects for the optical detection thermal emission 
from the "very hot Jupiters" in some detail and find the ob- 
servations promising, even from the ground. They calculate 
exposure times necessar y for detection using blackbody planet 
models and also those of iHubenv et al.l d2003l) . They highlight 
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OGLE-TR-56b and OGLE-TR-132b, the hottest planets yet de- 
tected, as the most promising candidates for secondary eclipse 
detections. 

6. PLANETARY CLASSIFICATION 

Previously ISudarskv et alj d2000i 120031) examined the atmo- 
spheres and spectra of giant planets as a function of orbital dis- 
tance, and classified these atmospheres in terms of the conden- 
sates that may be present. The pL Class that we suggest here 
overlaps the Sudarsky et al. Class IV and V planets, which are 
both dominated by alkali absorption. Class V planets are dis- 
tinguished by having silicate and iron clouds in the visible at- 
mosphere. However, it is not yet clear what effect these clouds 
will have on hot Jupiter spectra and energy balance, because 
particle sizes, cloud vertical distribution, and even composition 
(MgSiC<3 vs. Mg2SiC<4, which have very different optical prop- 
erties) are not known. Given that the derived P-T profiles and 
cloud condensation curves for pL Class atmospheres are some- 
what parallel (see Figure 2) there may be a relatively narrow 
region of irradiation level (and T e s) where clouds are present in 
the visible atmosphere, rather than down at hundreds of bars. 
The pM Class are planets with atmospheres hotter than those 
considered by Sudarsky et al. 

Although planetary classification based on dominant conden- 
sates is certainly attractive, in practice this may by difficult to 
implement. Clouds are often gray scatters over broad wave- 
length ranges, so it may be challenging to recognize their pres- 
ence and definitively identify the composition of a given cloud 
deck if one is suspected, particularly with data extending over 
a limited spectral range. In our own solar system the visible 
clouds of Jupiter and Saturn are surely dominated by ammonia 
ice, as predicted by equilibrium chemistry condensation, but 
there is as yet no spectroscopic confirmation of this widely held 
belief. 

A planetary classification scheme based on atomic and 
molecular absorption features may be more practical to imple- 
ment, and is in keeping with the tradition of basing classifica- 
tion directly on spectral features. This practice would be anal- 
ogous to the spectral MLT(Y) classification of cool dwarf stars 
and brown dwarfs. The T sp ectral type begins w hen methane is 
detected in the near infrared ( lKirkpatrickl2005l) . and it has been 
suggested that the arrival of ammonia absor ption in the near IR 
should delineate the arrival of the Y dwarfs ( Kirkpatr ickl2005h . 
Young, hot substellar objects and candidate planets directly im- 
aged at wider separations from their parent stars have and cer- 
tainly will co ntinue to be classifi ed in the standard MLT(Y) 
scheme (e.g.. IChauvin et al 1 120051) . Thus, extending this clas- 
sification scheme to all planets, near and far from their par- 
ents stars, is a reasonable extension of the status quo. The 
letters "M" and "L" are used to designate dwarfs at the stel- 
lar/substellar boundary according to their spectral features, and 
for the planetary objects with similar spectroscopic diagnostics 
as M and L dwarfs, we propose to use the spectroscopic classi- 
fication letters pM and pL. Classes pT and pY would be natural 
extensions to cooler temperatures. 

Planetary classification will surely be complex, as some plan- 
ets will be dominated by absorbed stellar flux while others by 
their own thermal emission, which will yield different P-T pro- 
files at a given T e ff. The wavelength ranges at which planets 
will be most easily observed will also be a function of stel- 
lar glare and the sensitivity of future observatories. Planets 
that are eccentric will straddle two ore more classes, as sug- 



gested in Figure 1 fo r HD 147506b and HD 17156b (see also 
Sud arsky et alj [20051 for predictions at larger orbital separa- 
tions). All of these close-in planets will have, at a minimum of 
complexity, a day/night temperature difference that will make 
planetary classification difficult. Obviously mass will be impor- 
tant as well since with decreasing mass one eve ntually leaves 
the gas giant regime. Hansen & Barman (2007) have pointed 
out an interesting correlation of the known transiting planets 
with Safronov number; they appear to fall into two distinct 
classes. This indicates that classification of planets surely will 
be complex and multi-dimensional. P erhaps we will eventu- 
ally find the lHansen & Barmanl ([2007) classes possess differ- 
ent atmospheric abundances due to different accretion histories. 
Thus we do not propose a detailed scheme, but rather suggest 
that the newly available transit data support a method of atmo- 
spheric classification based on molecular and atomic absorp- 
tion, rather than the presence or absence of particular conden- 
sates. 

7. DISCUSSION 
7.1. TiO Chemistry Revisited 

Of course the boundary between the proposed pM and pL 
classes could be somewhat indistinct, as the abundance of 
TiO/VO does not instantaneously drop to zero upon Ti/V con- 
densation. For TiO, Figure[2] shows that the mixing ratio falls 
by a factor of ten ^40-180 K after initial Ti condensation at 
millibar to bar pressures (Lee & Lodders, in prep). The spectra 
of early L dwarfs show a gradual weakenin g the TiO and VO 
bands over progressing cooler sub-classes dKirkpatrick et al.l 
1999; Kirkpatrick 2005); this points to transition planets with 
depleted but still significant levels of TiO/VO whose atmo- 
spheric temperatures are ^100-200 K cooler than the curve of 
Ti condensation at solar metallicity. HD 209458b is probably 
near the warm end of this transition region since HD 189733b 
shows no evidence for an inversion. We note that at these same 
temperatures, FeH and CrH are al so prominent abs orbers at 
wavelengths blue-ward of 1.3 fim dKirkpatrick! 12005), before 
being lost to Fe and Cr solid. Whether these molecules will 
eventually help in pM and pL sub-typing will depend on the 
quality of data available at these wavelengths. 

Since HD 209458b appears to have a temperature inversion 
even at temperatures where the abundances of TiO is expected 
to be waning significantly, it is natural to inquire whether the 
abundance of TiO, may differ from the prediction from chemi- 
cal equilibrium. TiO is typically the major Ti-bearing gas at the 
P-T conditions where Ti-bearing condensates are expected (see 
LodderillOOl. However, Ti0 2 gas is next in abundance, and 
other Ti-bearing gases are also present. With decreasing tem- 
perature (but before Ti-condensation), the TiO gas abundance 
decreases because abundances of other molecular Ti-bearing 
gases increase. The abundances of all Ti-bearing gases drop 
steeply when Ti-bearing condensates form. However, if con- 
densation is kinetically inhibited, the TiO gas abundances will 
drop much more gradually with decreasing temperatures, be- 
cause the TiO is consumed only by other Ti-bearing gases, not 
condensates. 

The abundance of Ti02 gas is of interest because T1O2 gas 
is a likely precursor of the Ti02-building blocks in Ti-bearing 
condensates. The net reaction TiO + H2O = Ti02 + H2 is 
independent of total pressure (there are the same number of 
molecules on both sides of the reaction). Under the relevant 
conditions here, the oxidation of the TiO radical is fast. For ex- 
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ample, at 1750 K and 0.01 bar, oxidation of TiO to the equilib- 
rium T1O2 gas abundance takes less than a minute. This means 
that TiO abundances are unlikely to be influenced by any verti- 
cal mixing processes in planetary or brown dwarf atmospheres, 
and equilibrium TiO and Ti02 abundances apply. 

The formation of perovskite, CaTi03, the expected conden- 
sate at total pressures < 0.02 bar in a solar composition gas, 
requires reaction of Ti02 with CaO. Preliminary calculations 
show that the oxidation of monatomic Ca to CaO gas proceeds 
on similar timescales as TiO to Ti02 oxidation. With respect 
to chemical changes by updraft mixing and quenching it must 
be emphasized that metal gas phase reactions are very fast so 
that equilibrium abundances established at higher temperature 
regimes will quickly re-adjust to the appropriate low tempera- 
ture equilibria during upwelling. Although the inhibition of the 
formation of Ti-bearing condensates should be kept in mind, 
another promising avenue for further work at this boundary may 
be further investigation into V-condensation. As discussed in 
§3.2, here we assume that VO is depleted with TiO, although 
strict chemical equilibrium would allow abundant VO at tem- 
peratures 200 K cooler than needed for initial Ti-condenation 
lLoddersld2002l) . 

At this time it would be unwise to rule out other absorbers 
such as non-equilibrium gases or condensates dr iven by photo- 
chemist ry, which we do not consider here. As [Burrows et al.l 
d2007bl) correctly point out, it may be possible for the abun- 
dance of a photochemicall y-derived ab s orber to scale as a func- 
tion of irradiation level. Liang et al. (2004]) have found that 
methane-derived hazes (which are present in the stratospheres 
of our solar system's giants) would not be stable for a variety of 
reasons, mainly because methane is not abundant at these high 
temperatures and lo w pressures, an d haze particulate conden- 
sation is prevented. IVisscher et alj ([2006) also briefly inves- 
tigated the ph otochemistry of C- and O-dominated species in 
HD 209458b. Ma rlev et all d2007[) suggest that photochemistry 
involving multiple abundant elements, including sulfur, could 
be important for these planets. Other authors (e. g. Yelle 2004) 
have investigated the photochemistry of the low density upper 
atmosphere of these planets as well. Addition investigations in 
this area are surely needed. 

7.2. Other Issues 

The models presented here are for day-side average plane- 
tary atmospheric structures and spectra. For pL Class planets, 
redistribution of absorbed energy will likely lead to somewhat 
cooler atmospheres and lower fluxes than predicted here. pM 
Class planetary atmospheres should be very close to a "no re- 
distribution model," that is hottest at the substellar point and 
becomes cooler toward the planetary limb. A no redistribution 
day-side is more lum inous than a uniform day-side average (see 
Barman et al. 2005^ who computed planet-wide average, day- 
side average, and no redistribution models of HD 209458b and 
TrES-1). Therefore our models are guides for understanding, 
and deviations from the predictions for specific planets are ex- 
pected. The pM Class planets could potentially be modestly 
brighter than we have shown. 

To our knowledge, no models for the thermal evolution and 
contraction of pM Class planets have included TiO/VO opac- 
ity in their atmospheres. Instead it has been assumed that Ti 
and V have condensed into clouds. If it is found that WASP- 
lb, TrES-4b, and other pM Class planets have the atmospheric 
properties that we describe, this will mean that evolution mod- 



els of these planets will need to be recomputed with new model 
atmosphere grids. As shown in §4.3, the atmospheric extension 
due to the hot stratospheric temperatures and strong TiO/VO 
opacity at opt ical wavelengths cou ld lead to larger measured 
radii (see also Bur rows et alj|2007bh . 

8. CONCLUSIONS: TWO CLASSES OF HOT ATMOSPHERES 

Though ID radiative-convective equilibrium model atmo- 
spheres we have addressed the class of atmospheres, the pM 
Class, f or which TiO and V O are extremely strong visible ab- 
sorbers (Hube nv et al.l2003l) . This absorbed incident flux drives 
these planets to have hot (^2000+ K) stratospheres. There- 
fore, these planets will appear very bright in the mid-infrared, 
with brightness temperatures larger than their equilibrium tem- 
peratu res. This is the ca se for HD 149026b (Harrington et al.l 
120071) . and HD 209458b iKnutson et all (|2007a|). In addition, 
these planets will have large day/night temperature contrasts 
because radiative time constants at photospheric pressures are 
much shorter than reasonable advective timescales. The hottest 
point of the planet should be the substellar point, which absorbs 
the most flux, leading to perhaps negligible phase shift between 
the times of maximum measured thermal emission and when 
the day side is fully visible. T his appears to be the situ ation for 
planet v And b, observed by Harrin gton et alj d2006l) . Given 
that its irradiation level is intermediate between HD 209458b 
and HD 149026b we find that this planet is pM Class. Due to 
the fast radiative times, the day side of these planets may have 
P-T profiles that do not deviate much from radiative equilib- 
rium models. Although atmospheric dynamics will surely be 
vigorous, winds will be unable to advect gas before it cools 
to space. For these planets, high irradiation, the presence of 
gaseous TiO and VO, a hot stratosphere, the location of the 
hottest atmosphere at the substellar point, and a large day/night 
temperature contrast all go hand-in-hand. 

On the other hand, we have shown that in the pL Class, domi- 
nated by absorption by H2O, Na, and K, photospheric pressures 
and temperatures prevail such that advec tive timescales an d 
radiative timescales are similar (see also ISeager et alj 12005). 
Since atmospheric dynamics will be important for the redistri- 
bution of energy, the consequences for the structure and ther- 
mal emission of these atmospheres will be quite complex. The 
efficiency of energy redistribution will vary with planetary irra- 
diation level, surface gravity, and rotation rate. pL Class planets 
will have smaller day/night temperature contrasts and measur- 
able phase shifts in thermal emission light curves that will be 
wavelength-dependent. In addition, these pla nets may show 
variab ility in secondary eclipse depth (e.g. Ra uscher et al.l 
2007a). However, without a better understanding of the dynam- 
ics it is difficult to make detailed predictions at this time. Sec- 
ondary eclipse depths should range somewhere between values 
expected for a "full redistribution" model and inefficient redis- 
tribution. The published secondary eclipse data for pL Class 
planets Tr ES-1 and HD 189733b are all consistent w ith this 
prediction dFortnev et al.l2005tlFortnev & M arlev 2007). In ad- 
dition, the 8 jim lig ht curves for 51 P eg b and HD 209458b 
dCowan etal] 120071) and HD 189733b dKnutson et alJ l2007bl) 
are consistent with this prediction as well. Examination of 
Figu red] shows that HP 1 79949b is a pM Class planet, and in- 
deed ICowan et al.l d2007l) found the largest phase variation it 
their small sample for this planet, but the unknown orbital in- 
clination makes definitive conclusions difficult. 

Additional observational results will soon help to test the 
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models presented here. We find that transiting planets WASP- 
lb, TrES-4b, TrES-3b, OGLE-Tr-lOb, and TrES-2b will be 
in the pM Class, along with non-transiters v And b and HD 
179949. The low-irradiation boundary of this class is not yet 
clear, and planet HD 209458b shows that temperature inver- 
sions persist to irradiation levels where TiO/VO are ex pected to 
begin being lost to condensation dBurrows et al.l2007bl) . At still 
lower irradiation levels, the limit ed data for HD 18973 3b lead 
us to conclude it is pL Class dFortnev & Marlevll2007h . Sec- 
ondary eclipse data for XO-2b, HAT-P-lb, and WASP-2b will 
be important is determining how temperature inverstions (and 
the TiO/VO abundances) wane with irradiation level. Just as in 
dM and dL stars, the condensation of Ti and V is expected to be 
gradual process, so we fully expect transition objects between 
the distinct pM and pL class members. 

We will soon have additional information that will help shed 
light on the atmosphere of H D 20945 8bThe 8 ( im light curve 
for HD 209458b obtained by ICowan etail d2007l) shows little 
phase variation. However, if our theory connecting TiO/VO 
opacity and temperature inversions to large day/night contrasts 
is correct, we expect to see a large variation. Soon H. Knut- 
son and collaborators will obtain half-orbit light curves for HD 
209458b at 8 an d 24 /urn. The quality s hould be comparable to 
that obtained by iKnutson et al.1 d2007bl) for HD 189733b, and 
will put o ur theory to the test. HD 147506, with an eccentricity 
of 0.517 dBakosetal.ll2007al). and HD 17156, with a n eccen- 
tricity of 0.67 dFischer et aUteOOHlBarbieri et alj|2007l) . will be 
extremely interesting cases as the flux they receive varies by 
factors of 9 and 26, respectively, between periapse and apoapse. 
They should each spend part of their orbits as pL class and 
part as pM class. This makes predictions difficult, but large 
day/night temperatures differences at apoapse are likely. 

There has recently been considerable discussion on the rel- 
ative merits of multi-dimensional dynamical models and ID 



radiative-convective model atmospheres for these highly irradi- 
ated planets. Both kinds of studies provide interesting predic- 
tions. While it could be claimed that ID radiative-convective 
models are unrealistic because they "lack dynamics," they do 
include very detailed chemistry, vast opacity databases, and ad- 
vanced non-gray radiative transfer, which all dynamics models 
for these planets lack. Given the large diversity in predictions 
among the various dynamical models, which use a host of sim- 
plifications, the next step will be combining dynamics and ra- 
diative transfer, an idea which has been mentioned or advanced 
by a number of authors dSeager etal]|200llBarman et al.l ;2005: 
Fort nev et al.ll2006aT: iBurrows et al.ll2006c iDobbs-Dixon & Lin 
120071) . Eventually we will be able to give up ID /-type pa- 
rameters to treat the incident flux in a more realistic fashion for 
these exoti c atmospheres . We rec ently star ted working toward 
this goal in lFortnev etaD ([2006a, see also Dob bs-Dixon & Lira 
2007); work continues, and we believe it will have a promising 
future. We think that the predictions we have made here with a 
ID model will provide a framework for understanding the ob- 
servations to come. Additional observations for both pL and 
pM class planets, along with additional theoretical and model- 
ing efforts, should further clarify our understanding. 
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FIG. 1 . — Flux incident upon a collection of hot Jupiter planets. At left is incident flux as a function of planet mass, and at right as a function of planet surface 
gravity. In both figures the labeled dotted lines indicate the distance from the Sun that a planet would have to be to intercept this same flux. Diamonds indicate the 
transiting planets while triangles indicate non-transiting systems (with minimum masses plotted but unknown surface gravities). Red color indicates that Spitzer 
phase curve data are published, while blue color indicates there is no phase data. The error bars for HD 147506 (HAT-P-2b) and HD 17156 indicate the variation in 
incident flux that the planets receive over their eccentric orbits. Flux levels for pM Class and pL Class planets are shown, with the shaded region around ~0. 04-0. 05 
AU indicating the a possible transition region between the classes. "Hot Neptune" GJ 436b experiences less intense insolation and is off the bottom of this plot at 
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FIG. 2. — Model P-T profiles for planets withg=15 ms~ 2 and 7j nt =200 K at various distances (0.025 to 0.055 AU) from the Sun. This 7i nt value is roughly consistent 
with that expected for a 1 Mj planet with a radius of 1 .2 Rj. Condensation curves are dotted lines and the curve where CO and CH 4 have equal abundances is dashed. 
The 0.1X Ti-Cond curve shows where 90% of the Ti has condensed out. However, even then TiO is a major opacity source. For none of these profiles has TiO/VO 
been artificially removed. The kinks in the 0.035 AU profile are due to interpolation diffifulties, as the opacity drops significantly over a small temperature range. 
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FIG. 3. — Total abundance-weighted atmospheric opacity at P=l mbar and T=1700 (blue), 1600 (green), 1400 (red) K. The hottest temperature, 1700 K, is warmer 
than the temperature for Ti and V condensation, which begins to occur at 1670 K. The intermediate temperature, 1600 K, is 33 K cooler than temperature at which 
Ti is 90% condensed. The optical opacity drops due to the removal of gaseous TiO and VO into Ti- and V-bearing solid condensates. By 1400 K TiO and VO bands 
bands have given way to prominent water bands and alkali lines. A strong water vapor band at 1.4 fim is readily seen at all three of these temperatures, and indicates 
the relatively modest change in absorption by water vapor with a 300 K drop in temperature. 
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FIG. 4. — Model P-T profiles for planets at 0.028 AU from the Sun. Effects of different values 7j nt and g are shown. Thick solid profile is for g=15 m s~ 2 and 
7] nl =200 K. Lower and higher dash-dot profiles are for g=15 m s~ 2 with 7j nt values of 150 K and 250 K, respectively. Lower and higher thin solid profiles are for 
7i nt =200 K with g values of 25 m s~ 2 and 9 m s~ 2 , respectively. 
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FIG. 5. — Incident net "flux" (erg g s fircC , solid line, left ordinate) in five model layers for a cloud-free pL Class model with TiO/VO removed, g=15 m s , 
at 0.05 AU from the Sun. The dotted line (right ordinate) illustrates the running integrated flux, evaluated from short to long wavelengths. The layer integrated flux 
is read at the intersection of the dotted line and the right ordinate. Note the logarithmic scale on the left. At 0.45 mbar, although the absorption due to neutral atomic 
alkalis are important, more flux is absorbed by water vapor. Heating due to alkali absorption becomes relatively more important as pressure increases. By 420 mbar, 
there is no flux left in the alkali line cores, and by 4.1 bar, nearly all the incident stellar flux has been absorbed. (After Marley & Mckay, 1999.) 
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FIG. 6. — Similar to Figure|5] but for thermal emission, from 30 to 0.26 /im, from long wavelengths to short. Negative flux is emitted. The dotted line is again 
integrated flux, evaluated from long to short wavelengths. The negative value of a layer integrated flux in Figure[5]equals the integrated flux here. In addition, the 
scaled Planck function appropriate for each layer is shown as a dashed curve. The temperature and pressure of each layer is labeled. Cooling occurs mostly by way 
of water vapor, but also CO. As the atmosphere cools with altitude progressively longer wavelength water bands dominate the layer thermal emission. 
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FIG. 7. — Incident net "flux" (erg g _1 s -1 /im _l , solid line, left ordinate) in five model layers for a pM Class model with g=\5 m s~ 2 at 0.03 AU from the Sun. 
Only the left ordinate scale is the same as in Figure[5] This pM Class model possesses a hot stratosphere. At 0.45 and 4.4 mbar absorption of incident flux (heating) 
is due almost entirely to optical bands of TiO and VO. Water vapor also absorbs flux as in the pL Class model, but its effect is comparatively swamped. Also similar 
to the model shown in Figure[5] by 420 mbar most flux at optical wavelengths has already been absorbed. 
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FIG. 9. — Both figures show, as a function of wavelength, the atmospheric pressure that corresponds to a given brightness temperature. Note the differences in the 
y-axes. Top: Planet at 0.03 AU which has a hot stratosphere induced by absorption by TiO/VO. Bottom: Planet at 0.05 AU that lacks a temperature inversion. The 
right ordinate shows the corresponding radiative time constant at each major tick mark from the pressure axis. Note that this right axis is not linear. The labeled gray 
lines at right indicate an advective wind speed that would be necessary to give an advective time scale equal to the given radiative constant. (See text.) For instance, 
in the top panel, at log P=-2.0, the radiative time constant is 2300 s, and an advection time of 2300 s would require a wind speed of 31 km s . 
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FIG. 10. — Pressure-temperature (solid lines, bottom x-axis) and pressure-T ra j profiles (dotted lines, top x-axis) for the models at 0.03 AU (pM Class, thick lines) 
and 0.05 AU (pL Class, thin lines, with TiO/VO removed). TiO/VO has been removed when calculating the 0.05 AU profile. Comparing this profile to the one at 
0.05 AU in Figure|2] where TiO/VO were not removed, shows slight differences around 1 bar, where equilibrium chemistry predicts a local increase in the TiO/VO 
abundances. 
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FIG. 12. — Left: HD 209458b atmospheric P-T profiles for a model without TiO/VO opacity (blue) and two models with TiO/VO opacity (green and red). The 
value of 7j nl is 200 K. Right: Day-side low resolution model spectra for t hese profiles, along with observational data. The Knutson et al. i2007q) data are shown 
from 3.6 to 8 /im. At 24 /im the lower point is from|E)erning et al. 12005) and upper point is from D. Deming et al. ( in prep). Solid lines show thermal emission 
while dashed lines show scattered stellar flux. The MOST la upper limit from ~0. 40-0. 68 fim from Rowe et al. 12007) is shown as horizontal black line. 
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FIG. 13. — Model P-T profiles for five planets. HAT-P-lb is likely pL class, which is the assumption here. The other four planets are pM Class. 
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FIG. 14. — Low resolution planet-to-star flux ratios for transiting planets shown in Pigure fTSI pM Class planets, due to their high incident fluxes are warmer 
and thus brighter. Flux ratios may approach 10^* in the red optical. Weak water emission bands are seen in the mid-infrared for the pM Class planets. The bars 
with arrows indicate the wavelengths at which thermal emitted flux becomes 10 times greater than scattered incident flux, for the HAT-P-lb model (black) and the 
TrES-3b model (red). pL Class fluxes red-ward of ~0.7 pm are dominated by thermal emission, not reflected light. pM Class fluxes red-ward of ~0.4 fim are 
dominated by thermal emission. 



